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(57) ABSTRACT 

Method and apparatus for estimating signal-to-noise ratio 
(SNR) y of a composite input signal e(t) on a phase- 
modulated (e.g., BPSK) communications link. A first 
demodulator receives the composite input signal and a stable 
carrier signal and outputs an in-phase output signal; a second 
demodulator receives the composite input signal and a 
phase-shifted version of the carrier signal and outputs a 
quadrature-phase output signal; and phase error 0 £ (t) con- 
tained within the composite input signal e(t) is calculated 
from the outputs of the first and second demodulators. A 
time series of statistically independent phase error measure- 
ments 0 E (t]), 0 E (t 2 ), .... QzrfU is obtained from the 
composite input signal subtending a time interval At=t (: ~t l 
whose value is small enough such that y(t) and o(t) can be 
taken to be constant in At. A biased estimate y* for the 
signal-to-noise ratio (SNR) y of the composite input signal 
is calculated using maximum likelihood (ML) estimation 
techniques, and an unbiased estimate y for the signal-to- 
noise ratio (SNR) y of the composite input signal is deter- 
mined from the biased estimate y*, such as by use of a 
look-up table. 

12 Claims, 3 Drawing Sheets 
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REAL-TIME SIGNAL-TO-NOISE RATIO 
(SNR) ESTIMATION FOR BPSK AND QPSK 
MODULATION USING THE ACTIVE 
COMMUNICATIONS CHANNEL 

ORIGIN OF THE INVENTION 

The invention described herein was made by an employee 
of the United States Government and may be manufactured 
and used by or for the Government for Government pur- 
poses without the payment of any royalties thereon or 
therefore. 

TECHNICAL FIELD OF THE INVENTION 

This invention relates generally to communications sys- 
tems and, more particularly, to methods for estimating a 
signal-to-noise ratio (SNR) of digital signals, such as phase- 
shift keyed (PSK) signals such as binary phase-shift keyed 
(BPSK) and quadrature phase-shift keyed (QPSK) signals, 
without a separate pilot signal being used. 

BACKGROUND OF THE INVENTION 

The transfer of binary data by carrier phase shifts, com- 
monly known as binary phase shift keying (BPSK), is a 
well-known technique and an alternative to amplitude 
modulation or frequency modulation data transfer tech- 
niques. It is however common for these various techniques 
to be combined, for example, a sub-carrier modulated by 
phase-shift keying used as a signal for the amplitude modu- 
lation or the frequency modulation of a main carrier. 

Phase shift keying (“PSK”) modulation is frequently used 
to transmit digital data. PSK involves shifting the phase of 
the carrier according to the value of the digital data. For 
example, in binary PSK (“BPSK”) the “zeros” in the digital 
data may be represented by a 1 80 degree shift in the phase 
of the carrier, while the “ones” in the digital data may be 
represented by no phase shift. Other degrees of phase 
shifting may be used. Quadrature PSK (“QPSK”) involves 
phase shifts of 0, 90, 180 and 270 degrees. PSK typically is 
referred to as “MPSK” where the “M” represents the number 
of phases. (The term “M-ary” is also commonly used.) 

After a transmitter sends an MPSK signal over the 
selected transmission medium (e.g., telephone lines or radio 
frequency waves), a receiver detects the phase changes in 
the signal. In order to do this accurately, the receiver must 
extract the unmodulated frequency and phase (commonly 
referred to as the reference frequency and phase) of the 
carrier from the received signal. Traditionally, phase-locked 
loop (“PLL”) circuits have been used to acquire carrier 
phase. PLLs are relatively easy to implement with either 
analog or digital technology and, in general, are considered 
to have good “steady state” performance. Receivers incor- 
porating digital signal processors (DSPs) are well known. 

The assessment of the propagation conditions that prevail 
on a communications link, whether it is a point-to-point 
terrestrial link or an earth-space link, is of vital importance 
for the optimal operation of the communications link. The 
significant indicator of the communications quality of the 
link is the signal-to-noise ratio (SNR) y. If the value of this 
quantity goes below a certain given threshold due, for 
example, to atmospheric conditions such as rain, the bit error 
rate on the link becomes unacceptable for reliable commu- 
nications integrity. Instantaneous (or real-time) knowledge 
of the dynamic behavior of the SNR y is thus essential for 
the optimal implementation of procedures to mitigate the 
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further degradation of y. In many instances, a separate 
propagation receiver (or pilot signal) is employed to acquire 
an associated signal transmitted by a beacon on a commu- 
nications satellite. Measurement of the fading conditions of 
5 this signal level is then extrapolated to that of the operational 
communications link. 

Many methods have been advanced to estimate the SNR 
using an active modulated communications channels. For 
example, the output of the receiver matched filter can be 
to sampled, i.e., the voltage level V s for an output symbol, and 
the value is compared to a pair of a priori determined voltage 
levels ±a,a>0o.<vV^. The statistical frequency of occur- 
rence n F of values which fail to fall within the interval 
[-a, a] is calculated as well as the total number n T of 
15 samples considered. The ratio n^/n^is related to the symbol 
error probability by 
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— = 2G Jk~ 
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where E s is the symbol energy (i.e., energy per symbol), N 0 
25 twice the noise power at the output of the matched filter and 
Q( . . . ) is the well-known Q-function. This relationship is 
then solved for the SNR y=E 5 /N 0 by using the inverse 
function Q -1 ( . . . ). Although this method has several 
shortcomings, e.g., limited dynamic range and sensitivity to 
30 automatic gain control variations and inter-symbol interfer- 
ence, it suffers from an irreconcilable defect; the inversion of 
the function Q( . . . ), necessary to obtain the estimate y of 
y, is mathematically correct only if dQ/dy^O. In the event 
that dQ/dy— >0, the problem of determining y becomes 
35 ill-posed. That is, a small error in the estimate of n.pJn T leads 
to a large error in the estimate y. In fact, the ill-posedness of 
this problem is the major source of the lack of dynamic 
range. 

As is well-known, for a problem to be well-posed, in the 
40 sense of Hadamard, it must meet the following criteria; (1) 
for each set of data, there exists a solution; ( 2 ) the solution 
is unique; and (3) the solution is stable, i.e., depends 
continuously on the problem data. If a problem does not 
meet one or more of these criteria, the problem is considered 
45 to be ill-posed. (See, e.g., A. N. Tikhonov and V. Y. Arsenin, 
Solutions of Ill-Posed Problems (John Wiley & Sons, New 
York, 1977).) 

An entirely different approach can be imagined which 
avoids the ill-posedness of the technique described above, 
so and thus tends to be more robust in the presence of mea- 
surement errors inherent in the sampling of the matched 
filter output. In the bi- (or binary) phase-shift keying (BPSK) 
case, the sampled voltage V s is related to the bi-phase signal 
amplitude ±A and the corresponding in-phase noise compo- 
55 nent N c by 

V s =K(±A+N € ) 


where K is a constant proportionality coefficient incorporat- 
ing receiver gain factors, etc. Remembering that the goal 
60 here is to obtain an expression for the SNR y=E s /N 0 =A 2 / 
2<j n 2 where o A , 2 =(N c 2 ) for zero mean, white Gaussian noise, 
the method endeavors to obtain this ratio solely from the 
measured values of V s . Thus, to separate the noise term, one 
65 would try to form the average <V 5 > using the fact that <N,,)= 
0. However, the random bipolar nature of the signal ampli- 
tude A also yields a zero average, giving <V 5 )=0. The 



US 7,190.741 B1 


3 

technique that is then adopted in this approach is to form the 
absolute value IV J of each sample and then forming the 
ensemble average giving 


(\V s \) = m±A + N c )\)*KA 


so long as the condition A»O n prevails. Additionally, the 
sampled values V s are used to compute the variance 
<V/ )=K 2 (A 2 +a N 2 ) 

Hence, using the former expression to rid to the A 2 term 
to give 

therefore allowing one to write 


A 2 <|V,I) 2 

7 2 <t\ 2«V2>-<|V S |> 2 ) 


solely in terms of the voltage samples V v . 25 

A major drawback of the technique described above is the 
formation and use of the absolute value of the random 
quantity V s ; such an operation can drastically change the 
statistical characteristics of the random variable. 

Even though using the absolute value of the random 30 
quantity V s was done to get rid of the bipolar nature of the 
communications signal amplitude, as will become evident 
from an understanding of the techniques of the present 
invention, as presented hereinbelow, using the absolute 
value of the random quantity V s will be shown to be 35 
mathematically faulty and, moreover, needless. As taught by 
the present invention (described hereinbelow), one can, and 
in fact should, incorporate the bipolar nature of the signal 
amplitude into a rigorous statistical analysis; this character- 
istic of the signal being just as important as its other aspects. 40 
In addition, the aforementioned absolute value approach 
requires the use of a predetermined bit-stream format within 
the communications data composed of a series of l’s, thus 
necessitating a synchronization with, for example, a pre- 
amble within the modulation format. This results in a further 45 
complication of its implementation and will be shown to be 
unnecessary. 

It is an underlying purpose of the present invention to 
formulate the rigorous statistical basis for the correct esti- 
mation of BPSK signal SNR from what is known about its 50 
behavior at the input and output of the receiver demodulator. 
Instead of employing tacit and unwarranted assumptions 
concerning the nature of the communication signal for 
analytical simplification, a complete consistent statistical 
description of a BPSK signal will be provided to which the 55 
well-known techniques of maximum likelihood estimation 
theory can be applied. By employing, rather than neglecting, 
all the subtitles of the statistics describing the BPSK signal, 
an unbiased estimation procedure will be derived that makes 
simple use of its inherent phase characteristics at the 60 
demodulator. In the following description, a preliminary 
review of BPSK signal representation will be given which 
will lay the foundation for the statistical connection between 
Gaussian noise and SNR. Once an appropriate probabilistic 
description is obtained that establishes a rigorous contact 65 
between SNR and the measured phase error of the BPSK 
signal entering the receiver demodulator, the methods of 
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maximum likelihood estimation theory will be used to 
obtain analytical expressions for biased and unbiased esti- 
mates of SNR from easily measured phase errors. Finally, 
the straightforward modifications needed at the demodulator 
to implement the required phase measurements will be 
given. It should be noted that the resulting SNR estimation 
technique is also applicable for a quarterary phase-shift 
keying (QPSK) demodulator simply by applying it to one of 
the BPSK arms with appropriate modifications for the SNR 
expression. 

PRIOR ART PATENTS 

Before proceeding to a description of the present inven- 
tion, the following U.S. Patents are cited as useful examples 
of prior art relating to receivers, modulation schemes, 
demodulators and/or SNR estimation, and are incorporated 
in their entirety by reference herein. In some cases, com- 
ments elaborating on the shortcomings of the prior art and/or 
contrasting the prior art with features of the present inven- 
tion are made. 

U.S. Pat. No. 5,317,599 discloses a method and circuit for 
detecting the carrier-to-noise (CN) ratio of a quarterary 
phase-shift keying (QPSK) signal. A binary signal including 
in the QPSK signal is reproduced by a converting circuit 
including a QPSK demodulator (FIG. 2; element 4) and a 
difference converter. (The QPSK demodulator 4 synchro- 
nous-demodulates the voice subcarrier to provide coherent 
component signal I and orthogonal component signal Q.) 
The binary signal is then error-controlled by an error detect- 
ing circuit and the detection of a double-error. The occur- 
rence number of the double-error is counted during a pre- 
determined time period wmch is long enough to compare 
with the occurrence time when two successive double-errors 
occur. A microcomputer stores a conversion table corre- 
sponding to the theoretical relation between the CN ratio and 
the occurrence number. The microcomputer computes the 
CN ratio from the detected occurrence number. The theo- 
retical relation is formulated based on the assumption that 
noise distributes isotopically around the QPSK signal in the 
signal space in accordance with the normal distribution. 
Although the method of this patent is related to the present 
invention, it is limited to a particular modulation method 
(i.e., is not generic) and the estimation procedure is ill-posed 
(not well-posed). 

U.S. Pat. No. 5,337,251 discloses a method of detecting a 
useful signal affected by noise. A measurement is taken of 
the expected S/N ratio (SNR) of this signal over a time slice, 
a measurement of the estimated white noise alone is taken 
over another time slice without the useful signal, the mean 
energy of the noise and of the noise-affected signal is 
calculated, in each of their time slices, the theoretical 
detection threshold is calculated, the ratio of these two 
energies is calculated, and the ratio is compared with the 
calculated threshold, this threshold being greater than 1 
(ideal threshold). This patent describes detection of a signal 
within a noisy environment, and an estimation procedure is 
employed to accomplish this, but the solution is not well- 
posed. 

U.S. Pat. No. 5,721,694 discloses a method and system 
for filtering a noisy input comprised of an information signal 
and a noise signal. No prior time domain value information 
or frequency domain parameter information on the informa- 
tion signal or noise signal is available. A whitening filter 
renders the noise part of a noise signal into a substantially 
white noise signal. A Fourier transform system provides an 
amplitude spectrum output and a phase output. A multiplier 
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squares the amplitude spectrum output thus providing a 
power spectrum output. A zero mean sub-system generates 
a zero mean power spectrum output from the power spec- 
trum output. A square root sub-system squares the zero mean 
power spectrum output thus providing a resultant output. An 
inverse Fourier transform sub-system inversely transforms 
the resultant output to provide a time domain resultant 
output. An inverse whitening filter provides a substantially 
noise free output from the time domain resultant output. 
Although this method uses sampling and statistical methods 
to attempt to provide a noise-free signal for communica- 
tions, it requires a specific receiver structure and does not 
address SNR estimation for the general case. 

U.S. Pat. No. 5,764,707 discloses a method for signal 
demodulation, based upon the assessment of an error statis- 
tic which is independent of the content of the data stream. To 
determine the appropriate sampling point within a data 
period, multiple sampling points are utilized to generate a set 
of alternative decodings of the data stream. Associated with 
each of these alternative sets is an error statistic which is 
computed based upon a known characteristic of the trans- 
mitted stream. For example, in Differential QPSK, the 
encoded phase change is known to be a multiple of 90 
degrees. The error statistic would be a function of the 
difference between each sampled phase and the nearest 
multiple of 90 degrees. This error statistic provides a quan- 
tifiable measure of quality and noise margin associated with 
each sampling point. The decoded stream associated with 
the error statistic of highest quality is then selected as the 
proper decoding of the transmitted data stream. This inven- 
tion is particularly well suited for processing small data 
packets, such as ATM. Although this method uses an indirect 
measure of SNR, it is used for direct signal reception and, 
therefore is part of a specific receiver structure. 

U.S. Pat. No. 6,061,339 discloses a fixed wireless loop 
system having adaptive system capacity based on estimated 
signal-to-noise ratio. The disclosed method is for operating 
a fixed wireless system (FWS) having a radio base unit 
(RBU) that uses a code division multiple access (CDMA) 
airlink for communicating with a plurality of subscriber 
units (SUs). The method includes steps of (a) estimating a 
signal-to-noise ratio (SNR) of the FWS; (b) comparing the 
estimated SNR to a threshold SNR value; and (c) one of 
allowing service to another SU or not allowing service to 
another SU based on the result of the step of comparing. The 
step of estimating preferably uses a normalized SNR, and 
further employs the use of a null PN code detection tech- 
nique. Although this method uses an estimated SNR to 
accomplish the adaptive operation of a digital communica- 
tions system, the estimation is not well-posed. 

U.S. Pat. No. 6,201,954 discloses a method and system 
for providing an estimate of the signal strength of a received 
signal. The signal strength of a received signal in a com- 
munication system is estimated by the measured signal 
strength of a desired signal which is transmitted from a 
source device to a destination device. The energy of the 
desired signal and the energy of the noise are measured or 
computed. The measured signal strength can be computed 
by dividing the energy of the desired signal with the energy 
of the noise. The signal strength of the desired signal can be 
estimated as the maximum likelihood estimate of the actual 
signal strength. The signal strength of the desired signal can 
also be estimated based on the expected value of the energy 
of the desired signal. The signal strength of the desired 
signal can further be estimated based on the expected value 
of the energy of the desired signal and accounting for a 
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predetermined bias. Although this method addresses the 
SNR estimation, the statistical procedure used is not generi- 
cally applicable. 

U.S. Pat. No. 6,233,439 discloses signal-to-noise estinta- 
5 tion of forward link traffic channel for fast power control. 
The process of estimating Eb/No (bit-energy-to-noise-den- 
sity ) of the received signal generates two streams of power 
control bits that are transmitted to the base station. One of 
10 the power control bit streams is generated under the assump- 
tion that the frame rate has not changed. The other power 
control bit stream is generated under the assumption that the 
frame rate has changed. The base station, knowing whether 
or not the frame rate changed, chooses the appropriate 
15 stream to use in controlling the base station transmit power. 
Although this method addresses SNR estimation, it is lim- 
ited in terms of its applicability. 

U.S. Pat. No. 6,317,456 discloses methods of estimating 
the signal-to-noise ratio (SNR) of a signal in an orthogonal 
20 frequency division multiplexed (OFDM) system include the 
steps of representing the signal as a complex symbol and 
rotating the complex symbol in a complex plane so that the 
signal component of the function lies substantially on the 
25 real axis of the complex plane and the noise function lies 
substantially on the imaginary axis of the complex plane. By 
obtaining the means squared value of the complex symbol 
and the variance of the complex symbol and dividing them, 
an estimate for the SNR of the signal is obtained. FIG. 2 of 
30 this patent is a block diagram of a receiver in an OFDM 
system for demodulating a signal having a noise component 
so that the SNR of the signal can be estimated. Although this 
method is very relevant in that it addresses the SNR esti- 
mation problem, it uses a mean-square estimation procedure 
35 which is exemplary of an ill-posed problem, and the method 
also requires post-processing. 

U.S. Pat. No. 6,137,824 discloses a method for estimating 
signal and noise quality, and a receiver. The highest value of 
a correlation vector is compared either to the second highest 
40 value or to the mean value of the other values. If a specific 
spreading code is used in the reception of the signal, the 
comparison forms a ratio of the signal and disturbance, the 
ratio being applicable particularly to combining signals from 
the diversity branches of a diversity receiver, whereby the 
45 method makes it possible to combine the signals even prior 
to detection. If the reception employs a spreading code with 
which no signal received by the receiver has been transmit- 
ted, it is possible by comparing the correlation values to 
form a noise quality estimate, an estimate value higher of 
50 which meaning that the receiver is receiving a signal and not 
just noise. Although this method addresses the SNR estima- 
tion problem, it requires a specifically prepared receiver 
structure and is therefore not generic. The estimation pro- 
cedure is not well-posed, and the method also requires 
55 pre-processing. 

U.S. Pat. No. 6,128,494 discloses signal-to-noise and 
interface estimator for analog voice channel. The residual 
phase deviation of the supervisory audio tone (SAT) and/or 
signaling tone (ST) signal in an analog voice channel 
60 cellular communication system is used to estimate the 
underlying signal-to-noise ratio (SNR) and signal-to-inter- 
ference ratio (SIR). The phase deviation is estimated using 
a two-stage autoregressive filter/variance estimator, which 
may be implemented in firmware in a digital signal proces- 
65 sor. This method estimates SNR on an analog, rather than 
digital communications link. The estimation procedure is 
not well-posed, and requires pre-processing. 
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SUMMARY OF THE INVENTION 

It is a general object of the invention to provide an 
improved technique for estimating the real-time signal-to- 
noise ratio (SNR) estimate 7 of a communications link. 

It is a further object to ease the operational and financial 
burden to terminal operation, by estimating y of the com- 
munications link by measuring various quantities of the 
operational communications link itself, without necessitat- 
ing the use of a separate propagation link which is indepen- 
dent of the communications link. 

According to the invention, method and apparatus are 
provided for estimating signal-to-noise ratio (SNR) y of a 
composite input signal e(t) on a phase-modulated (e.g., 
BPSK) communications link comprises providing a first 
demodulator receiving the composite input signal and a 
stable carrier signal, and outputting an in-phase output 
signal; providing a second demodulator receiving the com- 
posite input signal and a phase-shifted version of the carrier 
signal, and outputting a quadrature-phase output signal; and 
phase error G^jt) contained within the composite input signal 
e(t) is calculated from the outputs of the first and second 
demodulators. 

Further according to the invention, a time series of 
statistically independent phase error measurements GeIL), 
0£(t 2 ), . . . , 0 £ (U) is obtained from the composite input 
signal subtending a time interval At=t i .-t 1 whose value is 
small enough such that y(t) and oft) can be taken to be 
constant in At. 

Further according to the invention, a biased estimate y* 
for the signal-to-noise ratio (SNR) y of the composite input 
signal is calculated using maximum likelihood (ML) esti- 
mation techniques, and an unbiased estimate y for the 
signal-to-noise ratio (SNR) y of the composite input signal 
is determined from the biased estimate y*, such as by use of 
a look-up table. 

The present innovation provides a proper mathematical 
basis for such SNR estimation. The basic significance of the 
procedures disclosed herein make them applicable for a 
wide range of situations involving the digital transmission of 
data, including; 

terrestrial and space communication; 

cellular and hand-held telephones; 

planetary communications; and 

optical (including fiber) communications. 

The estimation techniques of the present invention are 
very useful in that they use, as basic inputs, the very entity 
that is desired to be assessed— i.e., the digital data stream 
which conveys the communication of information. The use 
of actual communications link in the assessment of the 
prevailing SNR makes for not only a more direct and simpler 
measurement of the propagation conditions, but also 
includes the potential degradations that may be introduced 
by the communications equipment itself — such effects can- 
not be captured by using a separate propagation link which 
is independent of the communications link. 

The inclusion of all the complications of the statistical 
nature of the communications data stream results in a simple 
and straightforward ‘recipe’ to deconvolve the SNR char- 
acteristics; the simplicity inherent in the approach is the 
result of fully utilizing the unimpaired mathematical nature 
which lies at the basis of digital communications. 

The invention therefore provides real-time, in situ signal- 
to-noise ratio estimation for the assessment of operational 
communications links, using the active communications 
channel (link). 
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Advantages of the SNR estimation technique of the 
present invention include that: 

the technique is passive in that it does not require any 
specifically prepared communications systems, signals 
5 or structures, is not limited to a particular modulation 
method, does not require a specific receiver structure, 
and does not require either pre- or post-processing; 

the technique is generic in that is applicable to any digital 
communications system that used phase shift modula- 
10 tion; and 

the technique is well-posed in that the mathematical 
procedure used is correct in terms of the ‘stability’ of its 
statistical solution, which means that the solution is 
impervious to the ever-present small variations in mea- 
15 sured input data. 

Other objects, features and advantages will become 
apparent in light of the following description thereof. 

BRIEF DESCRIPTION OF THE DRAWINGS 

20 

FIG. 1 is a simplified block diagram of a typical BPSK 
demodulator of the prior art, showing the relationship 
between signal parameters (between inputs and outputs); 

25 FIG. 2 is a vector diagram showing the relationship of the 
various signal and noise quadrature components of the input 
of the BPSK demodulator of FIG. 1, according to the prior 
art; 

FIG. 3 is a characteristic curve relating the biased estimate 
3 q y*, calculated from phase error measurements 0^) accord- 
ing to equation (Eq.30), to the unbiased estimate y, accord- 
ing to the invention; and 

FIG. 4 is a simplified block diagram of an embodiment of 
a BPSK demodulator implementing the techniques of the 
35 present invention. 

DETAILED DESCRIPTION OF THE 
INVENTION 

40 The present invention employs a rigorous statistical basis 
for the correct estimation of communication link SNR of a 
BPSK, QPSK (and M-ary) phase-modulated digital signal 
from what is known about its statistical behavior at the 
output of the receiver demodulator. Many methods to 
45 accomplish this have been proposed and implemented in the 
past, but all of them are based on tacit and unwarranted 
assumptions, and are thus defective. However, the basic idea 
is well-founded — i.e., the signal at the output of a commu- 
nications demodulator has convolved within it the prevailing 
50 SNR characteristic of the link. The acquisition of the SNR 
characteristic is of the utmost importance to a communica- 
tions system that must remain reliable in adverse propaga- 
tion conditions. The present invention provides a correct and 
consistent mathematical basis for the proper statistical 
55 ‘deconvolution’ of the output of a demodulator to yield a 
measure of the SNR. The mathematical analysis is based on 
the well-known statistical analysis of signal and noise flow 
through a generic M-ary demodulator and the techniques of 
maximum likelihood estimation theory. By employing, 
60 rather than neglecting all the subtitles of the statistics 
describing the stochastic phase-modulated signal, unbiased 
estimation procedures will be derived that makes simply use 
of the inherent signal phase characteristics at the demodu- 
lator output. The use of such techniques will alleviate the 
65 need and expense for a separate propagation link to assess 
the propagation conditions prevailing on the communica- 
tions link. Furthermore, these techniques are applicable for 
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every situation involving the digital transmission of data 
over planetary and space communications links. 

In the description that follows, a number of equations are 
set forth, and are generally identified, for purposes of 
cross-referencing, with the standard notation (Eq.x). 

Preliminaries of Signal Representation 

FIG. 1 illustrates atypical, simple BPSK demodulator 100 
of the prior art, showing the relationship between signal 
parameters at its input and output. 

A mixer (e.g., multiplier) 102 has two inputs: an input 
receiving (e.g., from a receiver block, not shown) a com- 
posite signal e(t) having signal s,(t) and noise n(t) compo- 
nents, and an input receiving (e.g., from a phase-lock loop, 
not shown) a stable carrier signal 2 cos(co 0 t). The output of 
the mixer 102 contains sum and difference components of 
the composite input signal and the carrier signal. To reject 
(suppress, eliminate) the sum component, and preserve the 
difference component, the output of the mixer 102 is passed 
through a low pass (bandpass) filter 104, the output of which 
is the in-phase output A,.(t)+N, (t) of the demodulator 100. 

The phase-modulated signal s,(t) at the input of the 
demodulator is defined by 


Si (/) = A(l)cos(ai 0 / + 0,(f)), 


m = 


2n(i— 1) 
M ’ 


1,2; At = 2 


(Eq.l) 


where ( 0 o is the angular frequency of the carrier wave, A(t) 
is its time varying amplitude and 0,-(t) is its bi-phase state. In 
terms of quadratures, the noise-free BPSK signal is given by 


Si(t) = /tp)[cos(wof)cos(ft(f)) - sm(a!of)sin(0,(r))] = (Eq.2) 

A(t)cos(uJot)cos(Qj(t)) = A c (t)cos(coot), 

A c (/) = A(r)cos(0,-(f)) 


where 

A c (0^(/)cos(e,(l))=±A(0 (Eq.3) 

This demonstrates that the noise free signal only has an 
in-phase component A c (t) since the bi-phase states are 
0j(t)=O°, 0 2 (t)=18O° thus relegating the sine factors to zero. 

Similarly, the noise n(t) associated with the input is given 
by 

n(f)=W(/)cos(ov+0\<O) (Eq.4y 

where N(t) is the time varying amplitude of the noise and 
0,/t) is the associated phase. 

In terms of quadratures, the noise is given by 


n(t) = At(/)[cos(o>oOcos(0«(d) - sin(wo0sin(0/dt))] = (ELq.5) 

N c (t)cos(u)Qt) - N s {t)sm(coot) 


where 

M c ( t)^N( t)c os ( 6 A -( f ) ) (Eq.6) 

and 


10 

showing that, unlike the noise-free BPSK signal component, 
the noise is characterized by both in-phase and quadrature 
terms. The composite signal e(t) at the input of the demodu- 
lator is given by 


e(t) = j;(r) + n(i) = [A,.(r) + A^pficosOjo/) + W s (/)sin(<Jor) (Eq.8) 


where use has been made of equations (Eq.2) and (Eq.5). 
Writing the composite signal e(t) in terms of its in-phase and 
quadrature components and employing equation (Eq.8) 
yields 

15 e(r)sE c (t)cos(u) 0 t)-E s sia((Oot) (Eq.9) 

where 


E c (t)=A c (t)+N c (t) (Eq.10) 

20 

and 


(Eq.llJ 

Writing equation (Eq.9) in the standard form in which the 
25 signal and noise components are written in equations (Eq.l) 
and (Eq.4), finally gives 

e(t)=E cos((o o (+0.e(;)) (Eq.12) 


where 

30 




e E (l) = tan 1 


-i 

E C 1 


(Eq.l 3) 


35 

Thus, 

E c (t)=E(t)cos(8Al)) (Eq.14) 


EM=E(l)sm^t)) (Eq.15) 

FIG. 2 illustrates, geometrically, the relationship of the 
various signal and noise quadrature components of the input 
45 of the BPSK demodulator in the case of the bi-phase signal 
s 1 (t)=A c (t)=A(t)cos(0 1 (t))=+A(t). In this figure, it is easily 
seen how the addition of the noise n(t) to the signal Sj(t) 
yields a total composite signal e(t) with a phase of 0^/t) with 
respect to the in-phase axis. The phase 0 £ (t) becomes, in this 
50 bi-phase signaling case, the phase error of the received 
signal Sjft), which represents the input to the demodulator. 

Given the preceding development of the relationships of 
the various signal components that enter into the demodu- 
lation of a BPSK signal, the problem which is being 
55 addressed by the present invention is the determination of 
the associated signal-to-noise ratio 


60 


y(t) = 


a 2 (P 


where the noise variance a N 2 is given by 


65 


2 , , (N;) (N;) 

°«(r) = 




(Eq.7) 


US 7,190,741 B 1 


11 

As described hereinbelow, this problem can be placed on 
a mathematically tractable basis and its correct solution rests 
on the statistical connection between the easily measurable 
parameters, not of the demodulator output, but rather of the 
demodulator input, i.e., E(t) and/or S^t), and the signal-to- 5 
noise ratio y. The fact that the demodulator input is signifi- 
cant is that it contains both the in-phase and quadrature 
components of the noise. Additionally, the bi-phase nature of 
the signal modulation, i.e., ±A(t), is a major impediment in 
the application of the prior art methods mentioned above, to 
According to the invention, this aspect of the signal must be 
properly accounted for in a legitimate treatment of the 
problem. It is therefore useful to examine (next section) 
what is known about the statistics of the noise process. 

15 

Noise Statistics and the Signal-To-Noise Ratio 

The noise components N c (t) and N s (t) are random quan- 
tities whose statistics are usually considered to be governed 
by a Gaussian random process. Additionally, it is known that 
such quadrature components are statistically independent. 20 
Letting o A , 2 (t) be the variance of each of these noise com- 
ponents, one has for their joint probability density 


p(NM, = 


( 1 ) N?(t) + Nl(t) 

V27ro^(oJ eXP 2tr 2 N (t) 


(Eq.16) 25 


Using equations (Eq.10), (Eq.ll), (Eq.14) and (Eq.15), 
one has in terms of the composite signal parameters E(t) and 30 
e^t), 

N c U)=E(t)cos(e^t)yA c (t) (Eq.17) 


and 

NJt)=E(l)s\a(6 f (t)) 


35 

(Eq.18) 


Transforming the probability density function of N c and 
N s (dropping the time dependence for notational clarity) 
given by equation (Eq. 1 6) into one that is a function of E and 40 
0 E (t) involves the transformation 


P (E, e E ) = P(N C , N s )\ 


3(N C , N s ) 
d(E, e E ) 


(Eq.19) 

45 


which requires the Jacobian 


3N C dN c 
\d(N c ,N s )\_ He Jiy 
I d (£, e E ) I “ 8N S dN, 

He HH 


50 

(Eq.20) 


55 


where equations (Eq.17) and (Eq.18) are used to obtain the 
indicated result. Thus, using equations (Eq.l6)-(Eq.20), one 
has in terms of the composite signal parameters 

, 60 


P(E, %) = 


(4) 

exp 

Ur) 

expj 


f EcosOf - A, ) 2 + /rsirrfy- 


2c rl 


(Eq.21) 


1 - 2A C EcosOe 
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The signal-to-noise ratio y characterizing the composite 
signal input to the demodulator is given for the BPSK case 
as y(t)=A(t)/o A , 2 (t). Writing equation (Eq.21) in terms of y 
and remembering equation (Eq.3) gives 


p(E, e E | y, ttj, <r N ) = 


(Eq.12 


' E 
.2-ncr\ 


rW-? 

N> 2 


is 2 g,).exp 


E 2 - 2Ecr N Hy cosB, cos8 E 
2o-% 


where the density function p(E,0 E ) is now written as a 
conditional probability density p(E,0£lY,0,-,a w ) governing the 
values of E and 0 E conditioned on the values of y, 0„ and a N . 
This is done with the hope of being able to obtain an 
expression, using equation (Eq.22), connecting the measur- 
able values of E and/or 0 E at the output of the demodulator 
only to the associated value for y , thus allowing one to 
statistically estimate y from such easily measured values. 

Having secured the above relationship governing the 
probability density of the envelope and phase of the com- 
posite signal, i.e., E and 6 E , and the prevailing values of the 
signal-to-noise ratio y, the variance of the noise power a N 2 , 
and the phase state of the transmitted signal 0,, one can 
obtain two relationships involving measurable statistics of 
E, or 0 E , and y. For example, one can obtain the statistics, 
conditioned on the value for y, that describe the composite 
signal envelope over all possible signal phase states simply 
by summing equation (Eq.22) over the two possible values 
for 0 ; (t) and integrating over all possible values of the phase 
0 E , using the following integral taken from I. S. Gradshteyn 
and I. M. Ryzhik, Table of Integrals, Series, and Products 
(Academic Press, New York, 1980). Eq.(3.534.2) (“Gradsh- 
teyn”): 


( E ) 

, 7 V ( E 2 | (Eq^) 




£ 

l-L 

(ircrj, 


VyE 

cru 

l“p(-|)exp| 


COS#£ 


+ exp 


VyE 

(Tn 


COS $£ 


d0 E = 


*)r 2co J^£ 


2cr\ jjo 
g) e xp(-|) e xp(-^)/„| 


V CTN 


cos 6e \d0E z 


£ 2 iJVrE _ 

ar N 


p(E | y, cr w ) 


which is a form of the well-known generalized Rayleigh 
distribution for the bi-phase signal envelope. Although this 
expression is useful for many purposes, it is not sufficient for 
the estimation of the quantity y from measurable statistics 
involving E(t) since, as can be seen from equation (Eq.23), 
it explicitly involves the parameter <s N 2 . Thus, using such an 
approach, one must know, a priori, the prevailing values of 
a N 2 . Further consideration shows that 1 his must be the case 
since the noise n(t) has two degrees of freedom which 
contributes to the composite signal envelope, i.e., N(t) and 
0,v(t), as seen from FIG. 2. Integrating over all possible 
values of the phase characterizing the demodulator output 
leaves the phase of the noise to be determined by knowledge 
of o^. Thus, if one desires to obtain an estimate of y from 
measurements of E or some statistic related to E, using the 
technique described above, one would also need an a priori 
estimate of a N , thus rendering useless a straightforward 
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estimate of y from easily measurable signal parameters. (In 
contrast thereto, it is the express object of the present 
invention to obtain an estimate of y from such easily 
measurable signal parameters.) 

However, if one considers statistics that describe the 
composite signal phase error 0 £ (l) by summing equation 
Eq.(22) over the two possible values for 0,(t) and integrating 
over all possible values of the signal envelope E(t), one 
obtains (e.g., using Gradshteyn, Eq.(3. 562.4)) 


this expression are still functions of time, 0 £ =0 £ (t). etc. 
Assuming that a time series of statistically independent 
phase measurements 0 £ (tj), © £ (t 2 ), .... ©jAD can be 
obtained from the composite signal input to the demodulator 
5 that subtend a time interval At=t £ — 1 2 whose value is small 
enough such that y(t) and o(t) can be taken to be constant in 
At, one can use equation (Eq.25) to obtain a probability 
density function conditioned on the series of phase mea- 
surements given by 
to 


ft p(£,fe|y,ft,o»d£ = 
1 1 


I Vr£ 


■'N 


^fy E \ . , , _ 

expl + cos 6 e + exp cos 6 e 

1 < r N J [ oat 


dE = 


E 1 ' ^ ■y E 

£exp| r- cosh cos Of \dE ~ 

1 2 (t 2 n ) <r N 1 


|ijexp(-|) J°° rexp|-yJcosh(V7rcos0 £ )<fr 


V7 cos(fe)exp(-7 s in 2 fe)erf| 


cos 6e 


(Eq.24) 


^r ex p(“?) = I r) 


p(y 1 fete), fete), - , fete)) = PI P ( r 1 fed;)) 

>=i 


(Eq.26) 
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by 


From this, one can form the likelihood functional defined 


20 


Uy I fete), fete), , fete)) = (Eq.27) 

k 

In(p(y Ifete), fete), ■■■ , fete))) = In (p(ylfete))) 

l=i 


25 


30 


where the third integral results from the change of variables 
mE/a N and erf( . . . ) in the fourth line is the ‘error function’ 
or probability integral defined by 


According to the method of maximum likelihood, the 
corresponding estimate y* of yis found from this functional 
as that value of y for which L( . . . ) is a maximum. Hence, 
the estimate y* of y which prevails over the set of phase 
measurements 0 £ (t,), o £ (t 2 ), . . . , 0 £ (t £ ), of a BPSK 
modulated signal is given by 


35 


dHy I fete), fete), "■ , fete)) 

fly 


(Eq.28) 


erf(x) s 


2 r 
y fit Jo 


exp(-t 2 )dr 


This formulation does not involve any a priori informa- 
tion other than the signal-to-noise ratio y. Thus, equation 
(Eq.24) gives a relationship involving the conditional prob- 
ability density of the phase of the demodulator output given 
a value for the signal-to-noise ratio y of the composite signal, 45 
with no other a priori information necessary. This expression 
is indeed a candidate for the basis of estimating y from 
knowledge of values for 0 £ . The proper mathematical foun- 
dation for this procedure will now be given. 

The Maximum Likelihood Estimation of Signal-to-Noise 50 
Ratio from Phase Measurements of the BPSK Demodulator 
Input 

Equation (Eq.24), giving the conditional probability den- 
sity governing the values of 0 £ given a prevailing value for 
y, can formally be written in the opposite sense, i.e., a 55 
conditional probability density governing the values of y 
given a prevailing value for 0 £ , viz., 

p(ylfe)= ™ 60 

Vy cos(fe)exp(-7sin 2 fe)rrf 

65 

Although not explicitly shown in equation (Eq.25) for 
notational simplicity, all signal components that enter into 


-cosfe +-exp(--) 



It is important to remember that the bipolar nature of the 
modulated signal as well as the usual assumptions of Gaus- 
sian noise are already convolved in the probabilistic descrip- 
tion of equation (Eq.25). Also, unlike the general case where 
y and a N 2 are functions of time, they can now be taken as 
constant and independent of time during sufficiently small 
measurement intervals At. A quantitative measure for ‘suf- 
ficiently small’ will be given below. 

It is now necessary to find the root y* of equation (Eq.28). 
Using equations (Eq.25)-(Eq.27) in equation (Eq.28), and, 
for analytical tractability, neglecting the last term of equa- 
tion (Eq.25) (an approximation which holds for large values 
of y), and performing the required differentiation yields 


fli 

dy 


sin 2 (0f(/i))++ 


1=1 

Yj ^{ ln j‘ ? ' / ( A /f COS(fe(,|)) |} 


(Eq.29) 


Completing the differentiation of the last term results in a 
rather unwieldy expression. Using the fact that the erf( . . . ) 
function tends to a constant for large values of the argument, 
i.e., for large y, one can neglect this term (consistent with the 
neglect of the second term of equation (Eq.25)) and find that 
in this case, the optimal estimate for y of a bi-phase signal 
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with additive Gaussian noise based on a series of k phase 
measurements 0 £ (t,) of the composite signal input is given 
by 


1 ^ 
1 


klj sin 


(Eq.30) 


, y » 1 


Because of the approximations involved with reglecting 
the second term of equation (Eq.25) and the third term of 
equation (Eq.29), this maximum likelihood estimate is 
biased toward large values of y. This bias must now be 
removed for small values of y. To this end, one uses the 
original expression of equation (Eq.25) retaining the second 
term since this will be appreciable for small y. However, to 
maintain analytical flexibility, the erf ( ... ) function of 
equation (Eq.25) will continued to be neglected. Finally, the 
analysis becomes more amenable if the inverse function 1/y* 
given by equation (Eq.30) is used. Given these consider- 
ations, one has (using Gradshteyn Eq.(3. 381.1) and use of 
the simplifying transformations embodied m Eqs.(8.356.1) 
and (8.359.4)) for the unbiased estimate y of y * over all 
possible values of the phase G^ft) 


TX; 


sm 2 (8 E )p(f\0 E )d8£ -- 


w-af 


sin 2 6 E dQ E ++ 



(Eq.31) 


30 


35 


Thus, the unbiased estimate y of y is related to the biased 40 
estimate y*, given by equation (Eq.30), through the non- 
linear relationship of equation (Eq.31). (It should be men- 
tioned that the ensemble average of sin 2 (0 £ ) is calculated 
over the phase interval -n/2SQ E Sn/2 since p('yl9 £ ) is peri- 
odic in 0,, with period jt. (Hence, phase values that differ by 45 
±it essentially correspond to the same SNR values and the 
averaging interval chosen in equation (Eq.31) implements 
the separate identification of these intervals.) Since there are 
no other parameters of the problem that enter into equation 
(Eq.31), this relationship is universal for BPSK modulation 50 
and is plotted in FIG. 3. FIG. 3 is a characteristic curve 
relating the biased estimate y*, calculated from phase error 
measurements 0 £ (t,-) according to equation Eq.(30), to the 
unbiased estimate y. From equation (Eq.31) one has the 
following limiting behavior: 55 


5-0 \y* / 2 7 

60 

Thus, for large y, the two estimates converge for the 
reasons explained above. In the case of small y, the biased 
estimate approaches 2 since in this case, the phase errors 
GHf,-) are randomly distributed within their range -jt/2S0 £ 65 
(t,-)§7t/2. Note that equation (Eq.31) does not monotonically 
approach the value of 0.5 as y— >0; it overshoots the 0.55 


level before it turns toward the 0.5 limit. This “defect” is due 
to the approximation made above in the neglect of the 
erf( . . . ) function in equation (Eq.29), and can be corrected 
if desired, based on the teachings set forth herein and 
additional analysis. 

The implementation of this process is straightforward. 
Considering now the composite signal input to the demodu- 
lator in terms of equations (Eq.9) and (Eq. 1 4)— (Eq. 1 5), 
instead of the previous representation of equation (Eq.8), 
and keeping in mind the goal of securing the values of the 
phase error 0 E (t), one arrives at the modifications shown in 
FIG. 4, which illustrates exemplary modifications (addi- 
tions) to the BPSK demodulator of FIG. 1 to obtain phase 
error in composite signal. 

FIG. 4 illustrates an embodiment of a the present inven- 
tion, implementing the techniques described hereinabove. 
The overall demodulator 200 includes two iterations 100a 
and 1006 of the demodulator 100 described with respect to 
FIG. 1. The one demodulator 100a has a mixer 102a with 
two inputs, receiving the composite signal e(t) and the stable 
carrier signal 2 cos(co 0 t) and outputting a signal, through a 
low pass (bandpass) filter 104a, which is the demodulated 
the in-phase output A c (t)+N c (t). 

The second demodulator 1006 also comprises a mixer 
1026 and a low pass filter 1046. The mixer (e.g., multiplier) 
1026 has two inputs: an input receiving (e.g., from the 
receiver block, not shown) the same composite signal e(t) 
which is input to the mixer 102a, and an input receiving a 
phase-shifted version of the stable carrier signal 2 cos(oo 0 t) 
which is input (not phase-shifted) to the other mixer 102a. 
The phase shift is effected by a functional block 106 which, 
in the case of BPSK, shifts the phase by 90 degrees. Other 
angles can be selected for other phase modulation schemes. 

The output of the mixer 1026 contains sum and difference 
components of the composite input signal and the carrier 
signal. To reject (suppress, eliminate) the sum component, 
and preserve the difference component, the output of the 
mixer 1026 is passed through a low pass filter 1046, the 
output of which is phase-shifted with respect to the output of 
the low pass filter 104a. 

In this example, the output (In-Phase Output) of the low 
pass filter 104a is 

£//)=£(Ocos(0^0)=A c (f)+tVc(O. (Eqs.lO, 14) 

and the output (Quadrature-Phase Output) of the low pass 
filter 1046 is 

£i/)=£(Osin(e^t))=Ay/). (Eqs.ll, 15) 

According to the invention, for purposes of estimating 
SNR of a BPSK signal, a QPSK demodulator structure is 
“evolved” from the original BPSK demodulator structure (of 
FIG. 1). In this manner, the quantities needed to obtain the 
phase error inherent in the demodulated signal can be 
captures. Thus it can be seen that one must form a QPSK 
demodulator from a BPSK demodulator so that, in addition 
to the In-Phase output of the BPSK demodulator, one now 
also has the associated quadrature-phase output. 

These two outputs of the low pass filters 104a and 1046 
are supplied to a functional block 108 which takes (calcu- 
lates) the arctangent (or alternatively the arc-cotangent) of 
the two signals (Es which is sine-based, and Ec which is 
cosine-based), the output of which is the phase error (noise) 
0 £ (t) contained within the composite input signal e(t), for 
sampling. This output of the functional block 108 is pro- 
vided to a functional block 110 which performs the calcu- 
lations set forth in equations (Eq.30) and (Eq.31). 
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As k sampled values of 0 £ (t) are obtained (how to 
determine a value for k is discussed below), one then 
implements equation (Eq.30) to obtain the biased estimate 
y* of the signal -to-noise ratio y. Once this has been obtained, 
equation (Eq.3 1 ) is employed and solved for the correspond- 5 
ing unbiased estimate y; this latter procedure can be effected 
by use of a ‘look-up’ table that represents the universal 
graphical behavior of FIG. 3. 

The functional blocks 108 and 110 are preferably imple- 
mented in software, and the mixers (102a, 102/;), low pass to 
filters (104a, 104/;) and phase-shifter (106) are preferably 
implemented in hardware. 

Determination of Phase Sampling Length to Assure Con- 
stancy of y and a N 2 During the Interval 

As mentioned above, y and On can be taken as constant 
and independent of time during sufficiently small measure- 
ment intervals At. It is imperative that the prevailing values 
for y and a A , 2 not evolve during the time period in which the 
k samples of composite signal phase 0 £ (t,) are formed to 2Q 
obtain the biased estimate y*. The fundamental constraint 
which must be addressed during the measurement process is 
the fact that y does not have spectral components that exceed 
20 Hz for an atmospheric channel. (This is typically the 
upper limit for atmospheric scintillation at communications 2J 
satellite frequencies.) Thus, if the data rate of the commu- 
nications link is R d bps, and if a measurement of the phase 
0 £ (t ( .) occurs for k consecutive bit intervals, one must satisfy 
the Nyquist sampling constraint 


This inequality serves to bound the number of phase 35 
samples. (Thus, the measurement interval At should prefer- 
ably be less than 14o second.) 

It will be appreciated by those with ordinary skill in the 
art that the methods of the present invention may be readily 
implemented using a computer processor and associated 40 
software. Thus, the circuitry of FIG. 4 may be integrated 
with a microprocessor associated with the respective 
receiver and having software to perform the above-described 
calculations. An application specific integrated circuit 
(ASIC) or chip set may also be constructed to perform the 45 
inventive methods, or a digital signal processor (DSP) may 
be designed to perform these methods. All of these embodi- 
ments and equivalents thereof are intended to be within the 
scope of the invention. 

As described hereinabove with respect to FIG. 4, a BPSK 50 
demod may be modified into a QPSK one so as to provide 
the quandature-phase output, from which both outputs are 
operated upon to obtain the phase error and finally, the 
unbiased estimate of the SNR. Based on the teachings set 
forth herein, one of ordinary skill in the art to which the 55 
present invention most nearly pertains may apply the tech- 
nique higher order phase modulated states. However, it is 
first important to note that the methods described herein- 
above apply only to a single class (albeit, a popular one, 
BPSK) of digital phase modulation schemes in which the 60 
digital information is encoded among M-possible evenly 
displaced phase shifts, each separated by 2*PI/M radians 
(see, e.g., equation (Eq.l)). For cases other than that of the 
BPSK case, one must go back into the theory and re-apply 
the analysis for the QPSK case, i.e., the four phase state case 65 
which requires equation (Eq.24) of the manuscript to be 
evaluated over four (i=l, 2, 3, 4) phase states. This will 
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result in a different functional form of the final expressions 
of equations (Eq.30) and (Eq.31). The same methodology 
would be applied to higher order (i.e., M>4) phase shifts. 
Thus, in all cases, the overall technique of the present 
invention and its morphology remain the same, but its 
application and form for each modulation scenario will be 
functionally (i.e., mathematically) different. 

Although the invention has been described in connection 
with various specific embodiments, those skilled in the art 
will appreciate that numerous adaptations and modifications 
may be made thereto without departing from the spirit and 
scope of the invention as set forth in the claims. 

For example, the simplest application of the techniques 
described hereinabove to QPSK and other M-ary phase 
modulated cases would involve using the method shown in 
FIG. 4 in conjunction with equation (Eq.30) with just one of 
the arms of the demodulator. In the interest of maintaining 
mathematical rigor, one could return to the analysis given 
above and employ, e.g., M=4 in equation (Eq.l) for the 
QPSK case. Thus, equation (Eq.24) would then employ a 
sum over i=4 possible phase states. It should be understood, 
however, that this would complicate the SNR estimate 
analogous to that of equation (Eq.30) for the BPSK case. 

Also, the invention is extendable to multi-level modula- 
tion techniques that include amplitude as well as phase 
modulation such as 16QAM (Quadrature Amplitude Modu- 
lation) which is a 4-level modulation technique transmitting 
four bits per symbol. 

Although the invention has been illustrated and described 
in detail in the drawings and foregoing description, the same 
is to be considered as illustrative and not restrictive in 
character — it being understood that only preferred embodi- 
ments have been shown and described, and that all changes 
and modifications that come within the spirit of the inven- 
tion are desired to be protected. Undoubtedly, many other 
“variations” on the “themes” set forth hereinabove will 
occur to one having ordinary skill in the art to which the 
present invention most nearly pertains, and such variations 
are intended to be within the scope of the invention, as 
disclosed herein 

What is claimed is: 

1. Method of estimating signal-to-noise ratio (SNR) y of 
a composite input signal e(t) on a communications link, 
comprising: 

providing a first demodulator receiving the composite 
input signal and a stable carrier signal, and outputting 
an in-phase output signal; 

providing a second demodulator receiving the composite 
input signal and a phase-shifted version of the stable 
carrier signal, and outputting a quadrature-phase output 
signal; 

calculating phase error ©^(t)) contained within the com- 
posite input signal e(t) from the outputs of the first and 
second demodulators; 

estimating a biased estimate y* for the signal-to-noise 
ratio (SNR) y of the composite input signal using 
maximum likelihood (ML) estimation techniques; and 

estimating an unbiased estimate y for the signal-to-noise 
ratio (SNR) y of the composite input signal from the 
biased estimate y*; 

taking a series of k phase measurements 0 £ (t,) of the 
composite signal input where k is a number of sampled 
values; and 

calculating a biased estimate y* for the signal-to-noise 
ratio (SNR) y of the composite input signal, over all 
possible values of the phase error 0 £ (t), according to 
the following equation: 
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^ siir(0 £ (/,)) 


y»l; and 

calculating an unbiased estimate y from the biased esti- 
mate y*, according to the following equation: 


| = J\in 1 (ti E )p(y I e E )de E ■■ 



sin 2 0£<i0£ + 


sin z (0f)cos(0£)exp — -siir(#£) \cL6e = 



2. A method, according to claim 1, wherein the carrier 
signal, as received by the second demodulator, is phase- 
shifted by 90 degrees with respect to the carrier signal 
received by the first modulator. 

3. A method, according to claim 1, wherein; 

the phase error 0 £ (t)) contained within the composite 
input signal e(t) is calculated by taking the arctangent 
(tan -1 ) of the outputs of the two demodulators. 

4. A method, according to claim 1, wherein the composite 
input signal is a bi-phase signal with additive Gaussian 
noise. 

5. A method, according to claim 1, further comprising: 

obtaining a time series of statistically independent phase 

error measurements 0 £ (t,), 0 £ (t 2 ), ■ • ■ , 0 £ (t*) f rom the 
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composite input signal that subtend a time interval 
At=t -t, whose value is small enough such that y(t) and 
a(t) can be taken to be constant in At. 

6. A method, according to claim 1, wherein the unbiased 
estimate y is determined from the biased estimate y* by use 
of a look-up table. 

7. Method, according to claim 1, further comprising: 
determining an unbiased estimate y from the biased 

estimate y* by use of a look-up table. 

8. Method according to claim 1, wherein: 

the ensemble average of sin 2 3 4 5 (0 £ ) is calculated over the 
phase interval -7t/2S0 £ Sjt/2. 

9. Method, according to claim 8, wherein the communi- 
cations link has a data rate R D bps, and further comprising: 

measuring the phase 0£<t,-) occurs for k consecutive bit 
intervals, and 

ensuring that following limitation is met: 
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Rd 

-j- > 40Hz. 


10. A method, according to claim 8, further comprising 
obtaining a time series of statistically independent phase 

measurements ©^(tj), 0£(t 2 ), . . . , from the 

composite input signal that subtend a time interval 
At=t^— 1 1 whose value is small enough such that y(t) and 
o(t) can be taken to be constant in At. 

11. A method, according to claim 1, wherein the commu- 
nications link is selected from the group consisting of BPSK, 
QPSK and other M-ary phase modifications. 

12. A method, according to claim 1, wherein the commu- 
nications link is selected from the group consisting of 
terrestrial and space communication, cellular and hand-held 
telephones, planetary communications; and optical commu- 
nications. 



